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I.  INTRODUCTION 


This  memorandum  describes  the  computational  techniques 
used  in  the  numerical  evaluation  of  a  solution  to  the  sound  pres¬ 
sure  field  resulting  from  a  plane  wave  incident  on  a  cylinder 
and  concentric  cylindrical  section..  The  analytical  development 
of  the  solution  is  given  in  a  separate  technical  memorandum  (1)  . 

In  the  numerical  evaluation  of  the  solution  to  the  baffle- 
transducer  sound  field,  values  of  Bessel  and  Neumann  functions  are 
required  for  a  large  range  of  order  and  argument.  Consequently, 
an  initial  requirement  for  the  numerical  evaluation  of  the  solu¬ 
tion  was  to  obtain  the  required  values  of  the  Bessel  and  Neumann 
functions.  Published  values  of  the  functions  are  only  available 
for  a  limited  range  of  order  and  argument.  Furthermore,  the  use 
of  the  asymptotic  expressions  available  for  evaluating  the  functions 
is  limited  to  small  ranges  of  order  and  argument  unless  a  large 
number  of  significant  figures  is  available.  Techniques  suitable 
for  evaluating  Bessel  and  Neumann  functions  on  a  small  digital 
computer  have  therefore  been  developed.  These  techniques,  which 
are  described  in  (2)  and  (3) ,  are  a  part  of  the  computational 
procedure  described  in  this  memorandum. 

For  the  purpose  of  numerical  evaluation  it  is  convenient  to 
divide  the  baffle-transducer  sound  field  solution  into  several 
quantities  which  may  be  computed  individually.  The  computations 
of  these  quantities  then  become  sub-routines  of  a  master  program 
for  evaluating  the  solution.  The  separation  of  the  solution  into 
these  quantities  is  described  in  Section  II  of  this  memorandum. 

The  coefficients  which  appear  in  the  solution  for  the  baffle- 
transducer  sound  field  arise  in  the  form  of  an  infinite  set  of 
inhomogeneous  equations .  For  the  numerical  evaluation  of  the  solu¬ 
tion  it  is  necessary  to  approximate  the  infinite  set  of  coefficients 
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by  a  truncated  set  of  coefficients.  This  procedure  is  discussed 
in  Section  III.  The  validity  of  the  approximation  is  supported 
by  demonstrating  the  convergence  of  the  individual  coefficients 
of  a  finite  set  to  stable  values  independent  of  the  size  of  the 
set  of  coefficients  evaluated  once  the  set  is  large  enough. 

In  the  last  section  some  numerical  results  for  the  solution 
are  presented.  Since  experimental  measurements  of  the  response 
of  scanning  sonar  equipments  to  an  incident  noise  field  are  given 
in  terms  of  self -noise  directivity  patterns,  the  sound  field  at 
the  transducer  computed  from  the  baffle-transducer  sound  field 
solution  is  presented  in  this  form. 

Appendix  A  contains  a  description  of  the  technique  used  in 
evaluating  a  finite  set  of  simultaneous  equations.  The  solution 
of  such  a  set  of  equations  occurs  in  the  numerical  evaluation  of 
the  solution  of  the  baffle- transducer  sound  field. 

Appendix  B  contains  a  description  of  the  computer  program 
used  for  the  numerical  evaluation  of  the  baffle- transducer  sound 
field . 

In  Appendix  C  a  description  is  given  of  the  technique  employed 
in  simulating  the  beam-forming  process  on  a  digital  computer. 
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II.  A  COMPUTATIONAL  FORIi  FOR  THE  SOLUTION 


In  this  section  the  form  of  the  analytical  solution  to  the 
baffle-transducer  sound  field  given  in  (1)  will  be  expressed  in 
a  form  suitable  for  computational  purposes.  The  data  to  be  ob¬ 
tained  in  a  numerical  evaluation  of  the  solution  ircludes  the  pres¬ 
sure  amplitude  and  phase  of  the  sound  field  at  various  points  in 
the  vicinity  of  the  baffle- transducer . 

The  form  of  the  solution  given  in  (1)  is 


Pi  *  y 


n=o 


j;(ka) 


N^(kr)]  cos  ne 


(1) 


where  pj^  is  the  sound  pressure  in  the  annular  region  between  the 
baffle  and  transducer  (see  Figure  1),  J^(kr)  and  N^(kr)  are  the 
Bessel  and  Neumann  functions  of  order  n  respectively,  the  primes 
indicate  differentiation  with  respect  to  the  argument,  k  is  the 
wave  number  of  the  incident  wave  (k  =  2n/X,  where  X  is  the  wave¬ 
length)  ,  and  is  determined  from  the  relation 


J^(ka) 

N^(ka) 


1^0 


j;(ka) 

Nj^(ka) 


j:(kb)  +  iN!(kb)  WJ 

j=0  J  J 


+ 


cp  ^  |€jiJjj(kb)  -  €jiJjJ(kb) 
j=0 


J^(kb)  +  iN^(kb)-> 
Jj(kb)  +  iNj(kb)J  U 


(2) 


This  equation  is  an  infinite  set  of  inhomogeneous  equations  of  the 
form 
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a._  and  are  constants. 

^n  n 


The  expression  for  the  pressure  given  in  eq.  (1)  is  a  complex 
quantity,  hence  the  magnitude  of  the  pressure  can  be  wricten  as 


cosnG] 


2 


1 


n 


(4) 


where  represents  the  real  part  of  eq.  (1)  and  is  the  imaginary 

part  of  eq.  (1).  The  phase  $  of  the  pressure  field  at  any  point 
can  be  written  as 


i 


tan 


-1 


^  cosn0 


Z 


cosnG 

n 


n 


(5) 


Ipj^l  and  §  are  the  quantities  that  are  to  be  evaluated. 


The  real  and  imaginary  parts  of  Pj^  (R  and  I)  are  determined 
from  the  expression  for  given  by  eq.  (2).  R  and  I  can  be 
written  as 


n 


K 


n 


n  n  n  n 


(6) 


and 


I_  =  K 


n 


n 


O  + 

n 


(7) 
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where 


K  = 
n 


PnJ„(kr)  -  a„N„(kr) 

a  ^  +  p  ^ 
n  '^n 


(8) 


=  J^(ka)  =  1/2  [J^_j^(ka)  - 


(9) 


Pn  =  N^Cka)  =  1/2  [N^_j^(ka)  -  N^^^(ka)] 


(10) 


and  are  determined  from  the  following  sets  of  simultaneous 
equations . 


V  0^  (A.X"  -  6'^B  p  )  -  Vd^  (ay"  -  6  Va  )  =  F„  (11) 

^  I  I  I  n  n'^n  I  I  I  nnn'  n 


I  /A  _  c 


Z®*  ■  <W  =  «n  (12) 


n 


{"} 


o  j  2  f  ^  ^  • « • 


Z  ^  Vt  {"}] 


t“l » 3  >  5 • • • 


(13) 
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1 
^  2 


Vt  {t} 


^=0,2,4,  . .  . 


Ill 

I  €,(-1)”  I- 

)3)5*  •  •  • 


t’’t  {"}  ] 


j 


A.  - 


2  .  Q  2 


4n^  [J^(kb)P^  -  a^N^(kb)] 

_  2  ^  2  ,  Q  2\ 


b. J. (kb)  -  a.N. (kb) 
T  —  J  J  J  J 

j  - - 

aj  +  bj 


a. J. (kb)  +  b.N. (kb) 
M  —  J  J  J  J 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
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a 


n 


b 


n 


=  j;(kb)  =  1/2  CJj^.i(kb)  -  Jn+i(kb)]. 

Sq  =  -  Jj^Ckb) 

(21) 

=  N;(kb)  =  1/2  [N^.i(kb)  -  N^^^(kb)]. 

b^  =  -  N^(kb) 

(22) 

«  |9  -f  sin.(n+j)£. 

n-j  p  Cn+j)cp 

(23) 

6^  =  1  ,  -t  =  n  and  ~  ^  ti  (24) 

The  computation  of  each  of  these  quantities  is  accomplished 
in  sub-routines  of  a  master  program  for  evaluating  the  solution 
given  in  eq.  (1).  The  master  program  and  sub-routines  are  de¬ 
scribed  in  Appendix  B. 
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III.  COMPUTATIONAL  TECHNIQUES 


The  primary  difficulty  in  the  numerical  evaluation  of  the 
baffle- transducer  sound  field  solution,  eq.  (1),  is  the  determination 
of  the  coefficients,  E^,  given  in  eq.  (3).  The  analytical  repre¬ 
sentation  of  is  an  infinite  set  of  inhomogeneous  equations  in 
which  each  coefficient  is  given  in  terms  of  itself  and  the  other 
coefficients.  A  qualitative  examination  of  the  functional  forms 
of  the  quantities  and  (see  equations  (2)  and  (3))  shows 

r  j  ^ 

that  the  factors 


{1} 


will  have  a  tendency  to  cause  the 


major  contribution  to  the  sum  on  i.  to  arise  for  values  of  i  in  the 
vicinity  of  n.  (These  quantities  have  convergence  characteristics 
similar  to  the  function  ^  •)  would  expect,  therefore,  that 
it  should  be  possible  to  obtain  accurate  numerical  values  for  the 
first  N  of  the  coefficients  E^  if  the  indices  in  the  summation  (n  and 
l)  in  eq.  (3)  are  restricted  to  a  finite  number  Ms  N.  Equation  (3) 


then  becomes  a  set  of  M  inhomogeneous  equations  in  M  unknowns  which 
can  be  evaluated  numerically  (see  Appendix  A) .  Equation  (1)  can 
then  be  evaluated.  The  number,  N,  of  accurate  coefficients  required 
is  determined  by  the  convergence  of  eq.  (1).  The  number,  M,  of 
coefficients  which  must  be  included  in  eq.  (3)  to  obtain  accurate 
values  for  the  first  N  of  the  E’s  is_ determined  by  the  convergence 
effects  of  the  factors 


{!}  •“  {:} 


This  technique  has  been  substantiated  by  determining  values 

of  the  coefficients  for  various  values  of  M,  assuming  a  fixed  value 

R  I 

for  N.  In  this  case  values  of  and  the  real  and  imaginary 

parts  of  E^,  were  computed  for  N  =  50  and  M  =*  100  and  M  =  150.  The 

results  of  this  comparison  are  given  in  Table  1.  It  is  seen  that 

I  R 

the  coefficients  and  have  reached  reasonably  stable  values 
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for  M  =  100,  i.e.,  for  M  =>  100  and  M  =  150  the  values  of  the 
coefficients  are  practically  the  same.  The  pressure,  eq.  (1),  was 
then  evaluated  using  each  of  the  sets  of  coefficients.  A  comparison 
of  the  values  obtained  is  given  in  Table  2.  It  can  be  seen  that, 
for  at  least  the  particular  set  of  parameters  used  in  this  evalua¬ 
tion  (k,  a,  b,  cp),  the  results  are  nearly  identical.  Hence,  the 
smaller  set  of  simultaneous  equations  (M  =  100)  is  satisfactory 
for  the  evaluation  of  eq.  (1).  The  pressure  amplitude  for  this  case 
is  shown  in  Figure  2. 

The  solution  was  verified  by  comparing  numerical  results  for 
the  case  of  cp  =  180°,  i.e.,  zero  baffle  width,  with  the  results  of 
an  evaluation  of  the  classical  solution  to  the  sound  field  resulting 
from  a  plane  wave  incident  on  an  infinite  length  rigid  cylinder. 

The  pressure  amplitudes  and  phases  obtained  from  the  two  solutions 
are  compared  in  Table  3*  In  addition,  the  pressure  amplitude  for 
the  case  of  zero  baffle  width  is  shown  in  Figure  3*  It  can  be  seen 
that  both  pressure  amplitude  and  phase  compare  quite  favorably  for 
the  two  solutions.  The  agreement  of  the  baffle- transducer  sound 
field  solution  with  the  simple  cylinder  solution  represents  a  sub¬ 
stantial  verification  of  the  techniques  used  in  the  baffle-transducer 
sound  field  computation. 
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TABLE  I 

A  COMPARISON  OF  THE  COEFFICIENTS  AND  COMPUTED  FROM 
TWO  DIFFERENT  SETS  OF  M  SIMULTANEOUS  EQUATIONS 

cp=165° 

ka=23.56  b=90" 


M=100  M=150 


On  oj 

and  dJ 
n  n 

Do 

.0350984 

.0422946 

.505417 

.492501 

-  .946866 

-.920642 

dJ 

.0566811 

.074510 

-.0635814 

-.080683 

D2 

-.898818 

- . 874366 

.810396 

.788585 

”3 

-.0689177 

-.0846728 

dJ 

.0998481 

.113703 

dJ 

.719360 

.700551 

D5 

-.615209 

-.599812 

D5 

.142809 

.154100 

D6 

-.202259 

-.209673 
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TABLE  I  (Cont'd) 


M=100  M=150 


and 

Dn  and  ®n 

4 

- . 508902 

-.496312 

“7 

.403891 

.393362 

-.271244 

-.275685 

^8 

.355985 

.356973 

.302783 

.293452 

-.223956 

-.214539 

.448126 

.446056 

®10 

- . 533861 

-  .529582 

-.168091 

- .  157640 

^11 

.178131 

.165982 

dJi 

- . 58O33O 

- . 575617 

4 

.657598 

.655479 

. 147571 

.137767 

^13 

-.159424 

- .  144676 

.642323 

.640644 
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TABLE  I  (Cont'd) 


M=100  M=150 


Dn  and  Dn 

Dn  and  Dn 

- .705866 

- .707561 

®14 

-.140504 

- . 126188 

.114500 

. 102767 

■>L 

-.685095 

-.689819 

°16 

.662697 

.667817 

°16 

.0750578 

.0671874 

.00130028 

.00585412 

.645939 

.650207 

°18 

-.419410 

- .418552 

“18 

.972592 

.973815 

Ul9 

-.174535 

-.177493 

°19 

-.726803 

-.725742 

^20 

.820809 

.817848 

d|o 

-.227778 

-.233440 

^21 

.239705 

.249350 

TRACOR,  INC. 


1701  Guadolup*  St 


Austin  I,  Tttxas 


TABLE  I  (Cont'd) 


M=100  M=150 


and  Dn 

Dn  and  dJ 

.927288 

.925269 

0^2 

-1.01166 

-1.01429 

D22 

.202999 

.216407 

“23 

- . 135361 

- . 148348 

D23 

-1.04896 

-1.05851 

^24 

1.13030 

1.14290 

.283495 

.282009 

.0245484 

.023998 

°25 

1.02915 

1.04665 

^26 

-.957032 

-.975116 

®26 

. 477974 

.474935 

^27 

-.105330 

-.0932969 

^27 

-.969128 

-.974884 

°28 

1.02553 

1.02458 

4b 

-.161015 

-.150995 
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TABLE  I  (Cont'd) 


M=100  M=150 


Dn  and 

Dn  and  Dn 

I)^ 

1)39 

. 164267 

.159353 

Dgg 

1 . 13682 

1.13368 

“30 

-1.10410 

-1.10283 

“30 

- .381412 

-.377503 

■'ll 

.0958926 

.100168 

■>^1 

-1.28357 

-1.28582 

U32 

1.18401 

1.18433 

0^2 

.276265 

.286122 

-.355106 

-.367154 

°33 

.962977 

.954309 

D34 

- .721675 

-.703262 

®34 

-.255979 

-.260411 

^35 

-.0542341 

-.0653916 

®35 

-.598181 

-.578178 

°36 

.291690 

222372 

“16 

-2.40690 

-2.42268 
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TABLE  I  (Cont'd) 


M=100  M=150 


Dn  and  dJ 

dS  and  dJ 

37 

.386058 

.409269 

°37 

1.09323 

1.10617 

D38 

-1.28874 

-1.31347 

“38 

.214077 

.221761 

D39 

.0858212 

. 0968643 

-1.33621 

-1.3577^ 

^40 

1.09312 

1.10185 

oJo 

.306737 

.326212 

D41 

-.217383 

-.238223 

“41 

.809808 

. 805783 

D42 

-.686704 

-.690052 

“42 

1.03805 

1.03123 

D43 

-.202972 

-.204207 

“43 

- . 947747 

-.915456 

D44 

1 . 13444 

1.10941 
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TABLE  I  (Cont'd) 


M=100  M=150 


P  T 

Dn  and  Dn 

Dn  and  D 

-.102192 

-.137711 

°45 

- .0708431 

-.0181885 

oJs 

1.12419 

1 . 14803 

°46 

-.983872 

-1.05459 

“46 

-.101903 

-.101561 

D47 

.0974928 

.189896 

“47 

-.831941 

- . 868575 

°48 

1.84956 

1.77326 

”48 

.767130 

.822620 

D49 

.133107 

.245294 

°49 

1.11337 

.989887 

Fig.  2  >  PRESSURE  AMPLITUDE  OF  BAFFLE-TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE 

REFERRED  TO  UNIT  AMPLITUDE  INCIDENT  WAVE  TRACOR,  D\WO  ^2000-178 

AUITIN,  TCIIM  •/(■/•I  WOM/ewe 


TRACOR,  INC. 


1701  Guadolup*  St 


Austin  If  Tsxas 


TABLE  II 

A  COMPARISON  OF  PRESSURE  AMPLITUDE  AND  PHASE  COMPUTED 
FROM  TWO  DIFFERENT  SETS  OF  SIMULTANEOUS  EQUATIONS 

cp  =  165® 

ka  =  23.56 
b  =  90” 


6 

m=i6o 

M=i56 

Amplitude 

Phase 

Pressure 

Phase 

3.75° 

. 18266 

315.073° 

.18429 

315.131° 

11.25® 

.45290 

140.602° 

.45404 

140.575° 

18.75° 

.47306 

318.353° 

.47510 

318.038° 

26.25° 

.36042 

130.972° 

.36218 

130.687° 

33.75° 

.25358 

272.848° 

.25616 

273.545° 

41.25° 

.27240 

61.582° 

.27226 

62.437° 

48.75° 

.35334 

248.045° 

.35656 

248.656° 

56.25° 

.54702 

95.670° 

.55156 

95.746° 

63.75® 

.89906 

287.212° 

.90664 

287.242° 

71.25° 

1.19226 

108.452° 

1.19968 

108.373° 

78.75° 

1.34960 

282.039° 

1.36010 

281.874° 

86.25° 

1.38304 

91.808° 

1.39450 

91.663° 

93.75® 

1.41814 

259.444° 

1.43102 

259.486° 

101.25° 

1.47368 

75.081° 

1.48828 

75.082° 

108.75° 

1.53762 

265.191° 

1 . 55680 

264.967° 

116.25° 

1.63716 

111.243° 

1.65322 

110.663° 

123.75° 

1.86178 

322.297° 

1.87308 

321.559° 

131.25° 

2.00686 

177-030° 

2.00810 

176.244° 

138.75° 

1.52218 

44.850° 

1.51294 

43.945° 

146.25° 

1.01020 

320.219° 

1.00212 

319 .754° 

153.75° 

1.10244 

257.328° 

1.09818 

256.083° 

161.25® 

1.39846 

225.442° 

1.36004 

224.134° 

168.75° 

1.47772 

178.390° 

1 . 44370 

177.096° 

176.25° 

1.22224 

210.162° 

1.19062 

208.776° 

Fig.  3  -  PRESSURE  AMPLITUDE  OF  BAFFLE-TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE  FOR  THE  CASE  OF  ZERO  BAFFLE 
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TABLE  III 

COMPARISON  OF  PRESSURE  AMPLITUDE  AND  PHASE  FOR  THE  CASE  OF 
ZERO  BAFFLE  WIDTH  OBTAINED  FROM  THE  BAFFLE -TRANSDUCER 
SOLUTION  AND  THE  CLASSICAL  SOLUTION  FOR  A  PLANE 
WAVE  INCIDENT  ON  A  RIGID  CYLINDER 


Angle  0 

Classical 

Solution 

Transducer-Baffle 

Solution 

Amplitude 

Phase 

Pressure 

Phase 

3.75° 

.0217562 

312.606° 

.0217520 

312.594 

11.25° 

.0665575 

131.703° 

.0665674 

131.702 

18.75° 

.115481 

309.974° 

.115477 

309.970 

26.25° 

.171593 

127.432° 

.171602 

127.426 

33.75° 

.238627 

304.140° 

.238633 

304.137 

41.25° 

.320976 

120.149° 

.320974 

120.144 

48.75° 

.423516 

295.508° 

.423511 

295.506 

56.25° 

.550975 

110.301° 

.550969 

110.297 

63.75° 

.706252 

284.682° 

.706250 

284.677 

71.25° 

.888185 

98.990° 

.888185 

98.983 

78.75° 

1.08940 

273.814° 

1.08940 

273 . 8O7 

86.25° 

1.29609 

90.058° 

1.29611 

90.049 

93.75° 

1.49085 

268.957° 

1.49085 

268.946 

101.25° 

1.65788 

92 . 060° 

1.65790 

92.048 

108.75° 

1.78789 

281.201° 

1.78788 

281.190 

116.25° 

1.87926 

118.4623° 

1.87925 

118.451 

123.75° 

1 . 93577 

326.111° 

1.93574 

326.100 

131.25° 

1.96302 

186.570° 

1.96303 

186.557 

138.75° 

1.97014 

62.340° 

1.97016 

62.328 

146.25° 

1.97546 

315.659° 

1.97549 

315.648 

153.75® 

1.99297 

227 . 878° 

1.99297 

227.868 

161.25° 

1.99995 

160.175° 

1.99996 

160.167 

168.75° 

1.99268 

114.654° 

1.99270 

114.649 

176.25° 

1.99787 

91.707° 

1.99790 

91.703 
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IV.  NUMERICAL  RESULTS 


The  numerical  evaluation  of  the  analytical  solution  to  the 
baff le-transducer  complex  has  been  completed  for  several  values 
of  the  parameters.  The  solution  Is  presently  being  studied  as  a 
function  of  Incident  wave  frequency  and  angular  extent  of  the 
baffle.  The  parametric  analysis  of  the  solution  will  also  Include 
the  effects  of  relative  baff le-transducer  spacing,  baffle  Impedance, 
and  transducer  diameter.  The  range  of  parameters  to  be  studied 
Is: 

1.  Incident  wave  frequency,  1  kc  -  15  kc 

2.  Angular  extent  of  the  baffle,  0°  -  120® 

3.  Baff  le-transducer  spacing,  1  ft.  -*  20  ft. 

4.  Baffle  surface  Impedance,  the  limiting  cases  of 
perfectly  rigid  and  perfectly  absorbing  surfaces  will 
be  analyzed.  If  experimentally  measured  values  of 
surface  Impedance  are  available,  these  can  also  be 
Included  In  the  parametric  analysis. 

5.  Transducer  diameter,  2.7  ft.  “*  l6  ft. 

The  response  of  a  sonar  to  the  anisotropic  sound  field  re¬ 
sulting  from  the  diffraction  of  a  plane  wave  by  the  baffle- transducer 
Is  obtained  by  (a)  computing  the  pressure  amplitude  and  phase  of 
the  diffracted  sound  field  from  eq.  (1)  at  points  on  the  transducer 
face  corresponding  to  transducer  elements,  and  (b)  using  this 
pressure  amplitude  and  phase  Information  as  Individual  transducer 
element  responses  In  a  simulated  sonar  beam  forming  process  (see 
Appendix  C).  The  beam  forming  process  has  been  made  an  optional 
sub-routine  In  the  master  program  for  evaluating  eq.  (1).  The 
numerical  results  of  an  evaluation  of  pj^,  eq.  (1)  are  then  obtainable 
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either  in  the  form  of  beam  responses  of  a  scanning  sonar  equip¬ 
ment  or  as  pressure  magnitudes  and  phases  at  points  on  the  transducer 
corresponding  to  element  faces.  A  parametric  analysis  of  the 
effect  of  the  baffle- transducer  diffracted  sound  field  on  sonar 
beam  forming  processes  is  then  easily  obtained. 

Some  numerical  results  for  the  total  diffracted  sound  field 
of  the  baffle- transducer  for  various  parameter  values  are  shovm 
in  Figures  4-23.  Both  the  pressure  amplitude  and  beam  pattern 
are  shown  for  each  set  of  parameter  values  studied.  For  these 
results,  the  transducer  diameter,  relative  spacing  of  the  transducer 
and  baffle,  and  baffle  surface  Impedance  are  each  constant.  The 
angular  extent  of  the  baffle  has  been  varied  from  0°  to  120®  (cp 
varied  from  180®  to  120°)  for  values  of  ka  =  8.77  (Figures  4  -  I3) 
and  17.54  (Figures  14  -  23). 
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Fig.  5  -  PRESSURE  AMPLITUDE  OF  BAFFLE-TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE 
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Fig.  6  -  PRESSURE  AMPLITUDE  OF  BAFFLE- TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE 
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Fig.  7  -  PRESSURE  AMPLITUDE  OF  BAFFLE-TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE 
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Fig.  8  -  PRESSURE  AMPLITUDE  OF  BAFFLE-TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE 
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Fig.  10 -SELF-NOISE  DIRECTIVITY  WVTTERN 
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Fig.  1 1  -SELF-NOISE  DIRECTIVITY  PATTERN 
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Fig.  12 -SELF-NOISE  DIRECTIVITY  RftTTERN 
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Fig.  14-  PRESSURE  AMPLITUDE  OF  BAFFLE-TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE 
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Fig.  15-  PRESSURE  AMPLITUDE  OF  BAFFLE-TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE 
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Fig.  16  -  PRESSURE  AMPLITUDE  OF  BAFFLE-TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE 
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Fig.  17  -  PRESSURE  AMPLITUDE  OF  BAFFLE-TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE 
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Flg.18-  PRESSURE  AMPLITUDE  OF  BAFFLE-TRANSDUCER  SOUND  FIELD 
AT  THE  TRANSDUCER  FACE 
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Fig.  19  -SELF- NOISE  DIRECTIVITY  WTTERN 
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Flg.20“SELF^  NOISE  DIRECTIVITY  BITTERN 
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Flg.2l  -SELF-NOISE  DIRECTIVITY  PATTERN 
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Flg.22 -SELF- NOISE  DIRECTIVITY  WTTERN 
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Flg.23 -SELF- NOISE  DIRECTIVITY  PATTERN 
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APPENDIX  A 

SOLUTION  OF  A  SET  OF  SIMULTANEOUS  EQUATIONS 

In  order  to  evaluate  the  solution  to  the  baffle- transducer 
sound  field  eq.  (1),  a  set  of  simultaneous  inhomogeneous  equations 
must  be  solved.  An  argument  has  been  presented  for  evaluating  a 
finite-set  of  equations  as  an  approximation  to  the  infinite  set 
presented  in  eq.  (3).  Here  it  is  assumed  that  a  finite  set  of 
simultaneous  inhomogeneous  equations  is  to  be  solved. 

The  finite  approximation  to  the  set  of  equations  given  in 
equations  (10)  and  (11)  can  be  represented  in  the  form 

^11^1  ^12^2  ^IN^N  " 

^21^1  ^22^2  *  ‘  *  ^2N^N  “  ^2 


^Nl^l  ^N2^2  '  ’ '  ^NN^N  “  ‘^N 
or 

N 

X  1  =  1,2.3- --N 

3-1 

Consider  any  particular  unknown  Xj.  Because  of  the  limited 
number  of  significant  digits  available  in  the  computer,  it  is 
desirable  to  choose  the  equation  which  depends  most  strongly  on  Xj 
to  eliminate  Xj  from  the  remaining  equations.  In  a  physical  sense 
this  corresponds  to  measuring  the  quantity  in  a  system  which  exhibits 
the  greatest  change  with  variations  of  a  particular  parameter  of 
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the  syatem  in  order  to  obtain  the  most  accurate  evaluation  of 
the  dependence  of  the  system  on  the  parameter.  Therefore,  the 
equation  in  which  the  magnitude  of  the  coefficient  of  Xj  is 
largest  as  compared  to  the  sum  of  the  magnitudes  of  the  rest  of 
the  coefficients  in  that  equation  is  chosen.  (In  actual  practice 
the  squares  of  the  coefficients  are  used.)  This  process  assures 
that  the  errors  incurred  in  the  additions  and  subtractions  neces¬ 
sary  to  the  evaluation  of  the  Xj's  will  be  minimized. 

In  order  to  solve  for  Xj  (initial  J  is  1) ,  the  following 
quantities  are  calculated  for  all  i's: 

N 

Pi  -  X 

A  value  of  i,  say  I,  is  then  chosen  such  that  the  quantity 

iful! 

Pi 


is  maximized.  After  I  has  been  determined,  the  quantities 


Ij 


^ij 


and 


Ct  =  — 

I  a 


IJ 


are  calculated  for  all  j.  Then,  for  ii4l  and  j  a  J,  the  quantities 

'ij 


®ij  “  ^ij 


Ci  =  c^ 


*1 

bi  a 
^iJ 


^IJ 

^  a 
au  IJ 


are  calculated,  where  a|j  =■  0,  and  ajj  =  a^j ,  j  <  J, 
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This  operation  reduces  the  coefficient  of  Xj  in  the  I  equation 
to  1  and  reduces  the  coefficient  of  Xj  in  all  of  the  other  equations 
to  zero.  Now  J  is  Incremented  by  1  and  the  process  is  repeated. 
After  this  process  has  been  carried  out  for  all  J  (l,2,3i«>*N) 
the  matrix  of  coefficients  has  one  and  only  one  coefficient  in 
each  column  equal  to  1  and  all  other  coefficients  in  that  column 
are  zero.  The  equations  are  then  re-ordered  along  with  the  new 
inhomogeneities  to  get  the  matrix  of  coefficients  in  the  form  of 
the  unit  matrix,  i.e.,  ones  down  the  diagonal  and  zeros  everywhere 
else  so  that, 

Z  aijXj  =  c; 

j=l 

reduces  to, 

a[iXi  =  c[  =  x^  , 

and  the  set  of  new  inhomogeneities,  c|, 
the  set  of  unknowns,  Xj^. 

We  were  originally  presented  with 
written  in  matrix  notation  as 

AX  =  C  , 

where  A  is  the  matrix  of  coefficients,  X  is  a  column  vector  of  the 
unknowns  and  C  is  a  column  vector  of  the  inhomogeneities.  We  have 
now  solved  for  X  and  obtained  a  solution,  X',  which  is  not  exact. 
Using  X'  and  the  original  set  of  coefficients.  A,  a  new  set  of 
Inhomogenelties  D  can  be  determined. 


becomes  the  solution  for 


a  problem  which  may  be 
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AX'  =  D  . 

This  equation  is  now  multiplied  by  a  constant  k,  which  is  to  be 
determined,  and  subtracted  from  the  original  equation  to  get: 

A(X  -  kX')  =  C  -  kD 

We  now  wish  to  pick  k  so  that  the  magnitude  of  the  vector,  C  -  kD, 
i.e.,  |C  -  kD|  is  minimized.  We  know  that, 

N 

|C  -kD|2  =  ^  (c^  -  kd^)2  =  y  , 

i=l 

and  to  get  y  a  minimum  with  respect  to  k,  y  must  satisfy 

or 

N 

^  2(c^  -  kd^)  (-d^)  =  0 

i=l 

N  N 

^Cidi-k  Y  =  » 

i=l  i=l 

and  it  is  seen  that. 


N 

I  "I'^i 


I 


i“l 
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Having  determined  k,  the  following  Is  calculated. 


x'^  =  c'^  =  kd  ; 

but 

*1  =  c| 
so  that 


The  x'^  are  taken  to  be  the  final  solution  to  the  set  of 
simultaneous  equations.  These  values  are  stored  on  magnetic  tape 
and  printed  out  along  with  the  c'^  and  c^.  The  c'^  and  c^^  are 
printed  out  in  order  to  determine  how  good  a  solution  the  x'j^  are. 
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APPENDIX  B 

PROGRAM  DESCRIPTION 

The  numerical  evaluation  of  the  baffle- transducer  sound 
field,  eq.  (1)  is  accomplished  using  a  master  control  program 
that  reads  the  required  sub-routines  into  the  computer  and  stores 
the  results  of  these  computations  for  use  in  subsequent  sub¬ 
routines.  The  evaluation  of  each  of  the  quantities  described  in 
Section  II  is  accomplished  with  a  sub-routine  of  the  master  con¬ 
trol  program.  A  flow  chart  for  the  computer  program  is  given  in 
Figure  IB.  A  description  of  the  various  sub-routines  is  given 
below.  Roman  numerals  identify  the  various  sub-routines  that 
appear  on  the  flow  chart. 

I.  Reads  k  and  x  (x=a  or  b)  from  paper  tape  and  calculates 
J^(kx),  N^(kx),  J^(kx),  and  N^(kx).  Stores  all  four 
sets  of  values  in  separate  files  on  magnetic  tape. 

II.  Reads  a^,  b^,  J^(kb),  and  N^(kb)  from  magnetic  tape 

and  calculates  L  and  .  Stores  both  sets  of  values 
n  n 

in  separate  files  on  magnetic  tape. 

III.  Reads  a^,  b^,  a  ,  and  from  magnetic  tape  and  calculates 
A^.  Stores  set  of  values  in  one  file  on  magnetic  tape. 

IV.  Reads  cp  from  paper  tape  r.i\d  calculates  •  Stores 

all  values  in  one  file  on  magnetic  tape.'" 

V.  Reads  L^,  b^,  a^,  and  {"}  from  magnetic  tape,  cp  from 
paper  tape  and  calculates  and  G^.  Stores  both  sets 
of  values  in  separate  files  on  magnetic  tape. 

VI.  Reads  J^(kb),  N^(kb),  a^,  and  from  magnetic  tape, 

cp  from  paper  tape  and  calculates  B^.  Stores  all  values 
in  one  file  on  magnetic  tape. 
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VII.  Reads  and  |^j  j  from  magnetic  tape  and 

calculates  x”  and  y”.  Stores  blocks  of  x”  and  y” 
alternately  in  one  file  on  magnetic  tape. 

VIII.  Reads  A  ,  B  ,  a  .  p  ,  x”,  y'^,  F  ,  and  G  from  magnetic 
Ti  n  n  n  ^  n  n  n 

tape  and  forms  matrix  of  coefficients  of  the  and 

n 

®n  columns.  Stores  each  column  as  a  block  on  magnetic 
tape  all  in  one  file.  Last  block  in  file  consists  of 
inhomogeneities  (F^  &  G^) . 

IX.  Reads  matrix  of  coefficients  from  magnetic  tape  and 

R  I 

solves  sec  of  simultaneous  equations  for  the  and  . 

R  I 

Stores  values  of  and  alternately  on  magnetic  tape 

in  one  file.  Also  prints  out  D^,  dJ,  F  ,  G_  and  F*  and 

n  n'  n'  n  n 

G^  calculated  from  the  solution.  (See  Appendix  A.) 

R  I 

X.  Reads  D^,  D^,  a^,  and  from  magnetic  tape,  k,  r,  and 

0  from  paper  tape,  and  calculates  |pj^l  and  i  for  desired 
values  of  r  and  0.  Prints  out  Ip^^I  and  $  and  if  directivity 
patterns  are  desired,  stores  |p^|  and  i  on  magnetic  tape 
in  one  file . 

XI.  Reads  (p^l  and  §  from  magnetic  tape  and  calculates  sum 
and  difference  self-noise  directivity  patterns.  (See 
Appendix  C.) 
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APPENDIX  C 

BEAM  FORMING  PROGRAM 


The  values  for  [p^^l  and  cp  which  are  obtained  from  eq.  (1) 
are  stored  on  magnetic  tape  so  that  they  may  be  used  as  input  data 
in  a  scanning  sonar  beam  forming  program.  The  values  of  Ipj^l  and 
cp  are  obtained  for  points  on  the  transducer  which  correspond  to 
transducer  elements.  The  beam  forming  program  duplicates  the  actual 
sonar  equipment  process  of  vectorially  summing  responses  of  the 
elements  to  the  right  and  to  the  left  of  the  center  of  the  beam. 
These  sums  are  then  added  and  subtracted  to  obtain  sum  and  dif¬ 
ference  signals.  The  magnitude  of  the  sum  and  difference  signals 
are  then  considered  the  beam  response  to  the  incident  sound  field. 

The  analytical  expression  for  the  sum  and  difference  signals 
is  written  as: 
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where  N  is  the  total  number  of  elements  used  in  forming  the  beam, 

P^  =  vectorial  sum  of  pressures  on  the  left  of  the  beam 

center, 

P^  =»  vectorial  sum  of  pressures  on  the  right  of  the  beam 
center, 

=  value  of  |pj^(  at  an  angle  0^  on  the  left, 

IPfl^  =  value  of  IPj^I  at  an  angle  0^^  on  the  right, 
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cp^  =  value  of  cp  at  an  angle  on  the  left 

D 

cp^  =  value  of  cp  at  an  angle  0^  on  the  right 

t  is  the  delay  time  for  an  element  at  an  angle  on  either 
side  of  the  center  of  the  beam  which  adjusts  the  phase  at  the  n’^ 
element  (at  the  angle  a  )  so  that  all  of  the  elements  appear  to  be 
in  a  straight  line  perpendicular  to  the  direction  of  the  beam. 


direction  of  beam 


a  cos  0  2na  cos  0 

uut^  =  2TTf  - - -  “  - -  =  ka  cos  0^  , 


where  c  is  the  speed  of  sound. 
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0  corresponds  to  the  direction  of  the  beam  lying  between 
the  elements.  The  number  of  elements  used  in  forming  the  beam 
may  be  varied  from  a  maximum  of  24  down  to  0 . 

In  addition  to  printing  out  the  values  of  IP^  ±  P^|  cor¬ 
responding  to  the  beam  being  directed  between  every  pair  of  elements 
on  the  transducer,  these  values  are  normalized  by  dividing  by 
the  number  of  elements,  N,  in  the  beam  and  these  values  are  printed 
out  in  db: 

|PL  t  db.  -  20  Xogio 

In  order  to  simplify  the  indexing  of  the  various  quantities 
involved  in  calculating  | ±  P^|  the  beam  is  taken  to  be  directed 
in  the  direction  0  =  90°  (See  Figure  1)  initially  and  successive 
values  correspond  to  0  being  incremented  by  J .5°  . 


TRACOR,  INC. 


I70l  Guadalupe  St 


Austin  I,  T«xas 


REFERENCES 


1.  TRACOR  Technical  Memorandum,  "An  Analytic  Solution  to  the 
Sound  Pressure  Field  Resulting  From  a  Plane  Wave  Incident 

on  a  Cylinder  and  Concentric  Cylindrical  Section,"  March  14,  1963 . 

2.  TRACOR  Technical  Memorandua,  "The  Evaluation  of  Bessel 
and  Neumann  Functions  Through  the  Use  of  Taylor  Series 
Solutions  to  Bessel's  Differential  Equation,"  December  3,  1962. 

3*  TRACOR  Technical  tfemorandum,  "Computation  of  Bessel  and 
Neumann  Functions,"  December  11,  1962. 


